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Abstract 

At the HERA collider at DESY, high energy electron and positron beams interact with 
proton beams. A review is presented of the variety of ways in which these collisions 
produce final states containing charm and beauty quarks. 



1. Introduction 



1.1. The HERA collider 



The HERA collider is the only experimental facility in the world where high-energy 
electron-proton and positron-proton interactions can be studied. Electrons (or 
positrons) and protons are accelerated in the DESY ring complex to energies of 27.5 
GeV and 920 GeV respectively (fig. |l]a) . The circulating beams are then brought to 
a focus at the two collision points where the HI and ZEUS experiments are located. 
In most of the e^p interactions, the incoming interacts with the proton by first 
radiating a virtual photon (fig. ^). The virtuality of this exchanged photon 
(for definition see section 1.5) lies within a wide range; its minimum value is limited 
by the kinematics of the radiation process to be of the order of the electron mass 
squared, while its maximum value is limited to 10^ GeV^ by the energies of the 
incoming beams. 

During the period 1992-2000, the HI and ZEUS experiments each accumulated 
approximately 135 pb~^ of integrated luminosity, which has enabled a broad range 
of physics to be investigated. At low , the exchanged photon resembles a real 
photon, and we shall usually refer to such collisions simply as 7p or "photopro- 
duction". Real photons, as will be discussed below, can display hadronic as well 
as electromagnetic properties. The hadronic behaviour falls away over a transition 
region where is of the order of 1 GeV^ , above which we are in the regime of deep 
inelastic scattering (DIS). Here the photon is primarily an electromagnetic object 
which couples to quark-antiquark {qq) pairs. At very high the exchange of Z and 
W bosons also becomes important. The 7*p luminosity is approximately uniform 
in logQ^, giving a prolific production of low-Q^ photons. 

Differences between e'^p and e~p collisions are of significance in some new par- 
ticle searches and at the very high values where electroweak effects occur. In 
the physics reviewed here, however, the sign of the incoming e beam is nearly al- 
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Fig. 1. (a) The HERA collider complex at DESY, showing the ep ring and the lower energy 
accelerator PETRA. The HI and ZEUS experiments are in the North and South halls respectively. 
The proton beam energy was raised from 820 to 920 GeV at the end of 1998. (b) ep interaction 
mechanism through the radiation of a virtual photon. The broad arrow denotes the final-state 
products. 

ways unimportant. For convenience, we shall normally refer to the incoming as a 
positron, including electrons by implication, since most of the integrated luminosity 
was taken with positrons. The symbol 7* will denote any virtual photon, while 7 
will refer, where appropriate, to the particular case of low-Q^, quasi-real photons. 

1.2. QCD physics with heavy quarks at HERA 

In the present survey, we consider those processes at HERA in which heavy flavours 
are produced, namely the production of particles containing charm (c) and beauty 
(b) quarks. HERA is par excellence a laboratory for the study of quantum chromo- 
dynamics (QCD), our present theory of the strong interaction. While e+e~ colliders 
such as LEP have enabled the production of qq pairs to be studied in depth, along 
with the associated radiation of gluons, the presence of the proton in the initial state 
at HERA allows a wide variety of further QCD processes to be investigated. The 
incoming 7*, with its range of virtualities, probes powerfully the partonic properties 
of the proton and the reactions initiated through the quarks and gluons of which 
the proton is a source. 

Heavy quark production can also, of course, be studed in hadron-hadron colliders 
such as the Tevatron. However a virtual photon is a cleaner probe than a second 
high-energy hadron, and it is experimentally possible to trigger on photon-induced 
processes that occur down to lower transverse momenta, down to a few GeV. 

The behaviour of QCD processes is governed by the running coupling constant 
Qfs, whose value is «0.12 at momentum transfers corresponding to the Z mass. 
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though significantly higher at a few GeV. Consequently, the lowest-order (LO) per- 
turbative calculations, which suffice to a first approximation, require frequently to 
be replaced by next-to-leading order (NLO) calculations or higher. This can present 
formidable technical difficulties, which can limit our experimental understanding of 
QCD. With light quark systems, a certain impasse may by now have been reached, 
but the consideration of 6 or c quarks presents a number of new factors. These arise 
because the production of a heavy object requires a corresponding momentum trans- 
fer which may be significantly higher than the QCD scale parameter Aqcd ~ 300 
MeV. Such processes offer a different perspective on the operation of QCD. The 
measured final-state hadron containing a heavy quark relates directly to the prop- 
erties of the QCD-governed hard scatter, and although the presence of a second 
momentum scale can complicate the theory, it provides a means of calculating cross 
sections using perturbative approaches that are unavailable, or less reliable, when 
dealing with exclusively light quark systems. 

There remain a number of difficulties, however. The effective mass of the c 
quark is uncertain within the range 1.3 to 1.5 GeV, a value which is itself not over- 
whelmingly higher than Aqcd- These facts translate into significant uncertainties 
in the perturbative calculations, comparable to the effects of higher-order QCD di- 
agrams. As will be seen, this makes the interpretation of many experiments less 
than straightforward. The higher mass 4.75 GeV) of the b quark should reduce 
these problems; nevertheless it will be seen that even the production of b hadrons 
currently presents challenges to theory. 

A major aspect of heavy quark physics concerns the weak decays of the hadrons 
that contain these quarks. Such decays may be studied at any collider that is able 
to produce the relevant hadrons in sufficiently large numbers. Many of the most 
interesting decay channels have low branching ratios, however, and HERA has not 
yet produced data in enough quantity to be competitive here; CP violating effects, 
especially, lie in the domain of the B factories and the Tevatron. The present survey 
therefore concentrates on the production mechanisms of heavy flavour systems, an 
area where HERA is able to play a prominent role. The author's aim is to provide an 
account which is accessible to non-specialist readers and gives an up-to-date picture 
of the state of this area of HERA physics. Experimental results are preliminary 
where indicated. 

1.3. The HI and ZEUS detectors 

The HI and ZEUS detectorsS'i are situated respectively at the north and south 
collision regions of HERA. Each detector comprises of a cylindrical central tracking 
chamber, surrounded by calorimeters. Both central tracking chambers are of drift 
chamber design. In the case of HI, the calorimeters consist of lead and steel con- 
verters in liquid argon, and are mainly located inside the solenoidal magnet which 
provides a magnetic field of 1.15 Tesla. In ZEUS, the calorimeter is a uranium- 
scintillator sandwich. It lies outside the solenoid, within which there is a field of 1.4 
Tesla. The return yokes are part of the muon detection system, with muon detec- 



Heavy Flavour Physics at HERA. . . 4 




P (d) 

Fig. 2. Examples of different types of photon-proton interaction: (a) direct (b) resolved (c) 
anomalous (d) elastic diffractive. 



tion arrays outside them. HI has a silicon vertex detector array; ZEUS is currently 
commissioning its own silicon system. 

Both detectors have a design which is approximately cylindrically symmetric 
around the beam line, but with a forward-backward asymmetry. In the forward 
(proton beam) direction there is a need for good hadron calorimetry in order to 
detect energetic particles associated with the proton, its remnants, and particles 
that are in general boosted in this direction. In the rear direction, the best possible 
electromagnetic calorimetry is required to measure accurately the deep inelastic 
scattered positrons. 

In both experiments the calorimetry is designed to be as hermetic as possible, 
with holes of a minimum size for the beam pipes. 

1.4. Types of photon-hadron interaction 

Since the ep processes of concern here are mediated almost entirely by photons, 
some basic classes of photon-proton interactions must be considered. The most 
important of these are illustrated in fig. ||, where the first three diagrams depict 
different hard QCD scattering processes. Direct processes (a) are those in which 
the entire photon couples to a quark-antiquark line at high transverse momentum 
Pt- Alternatively (b), the proton may interact through an intermediate hadronic 
state which is a source of partons that can scatter off those in the proton. These 
are known as resolved photon processes, and here the photon possesses a parton 
structure like any other hadron. The concept of a resolved photon can be extended 
to situations where the photon couples to a relatively low-p^ quark pair (c). These 
perturbatiyely calculable processes are sometimes referred to as anomalous photon 
processes.B 

Ultimately, the outgoing partons hadronise as jets or beam remnants. The sig- 
nature for a resolved or direct photon process at LO is the presence or absence. 
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respectively, of a photon remnant in the event. The direct and resolved classes of 
photon interaction are fully distinct only in LO processes such as those illustrated. 
At NLO the two classes can still be defined for calculational purposes, but a contin- 
uum of diagrams exists connecting them. In photoproduction, the resolved processes 
occur at higher rates than the direct process, but as the of the photon increases, 
the resolved cross sections decrease and the reactions become predominantly direct 
in nature. 

The fraction Xj of the incoming photon energy that enters a given QCD sub- 
process is an important parton-level quantity. In direct processes it is by definition 
unity, while in resolved processes it can take any value in the range (0,1). The 
parameter x^ is clearly defined in LO diagrams but at higher order there can be 
ambiguities. Experimentally the hadronic final state is measured, and a suitable 
final-state quantity is required which correlates with x-y in a chosen theoretical de- 
scription of the process. For dijet final states, two such experimental estimator^S 
are ^ 

„OBS _ ^jetsEJp e~'' ^ HjetsjE - Pz) 



XX — — ana x. 

^7 



^ 2E^ ^ ^,,entiE-pz)' 



Here as elsewhere Et denotes transverse energy, namely Esm9 for a single particle, 
summed over the particles in a jet. The laboratory pseudorapidity 77 is given by 

e 

77 = — in tan — , 

where 6 is the polar angle and the proton beam defines the forward {+z) direction. 
The proton has E~pz = to an excellent approximation, and a quasi-real incoming 
photon has E — pz = 2E^ — Y,event[E — Pz)- Snowmass — i.e. £^T-averaged — 
quantitiesB are used in x°^^ , while in x^™'^"" the sums are over all particles in the jets 
and in the event as a whole. One notes that for a massless object, Exe^^ = (E—pz); 
the two estimators differ in that x?^^ neglects jet-mass effects. The estimator x^^^ 
has been mostly used in HERA analyses. 

A further class of processes in which bb and cc systems are produced comprises 
dijfractive processes. An example is illustrated in fig. |2|i, where the photon forms an 
intermediate vector meson state V which scatters off the proton by the exchange of 
a pomeron. Both quasi-real and highly virtual photons may take part in diffractive 
processes at HERA; these will be the subject of Section 6 of this paper. 



1.5. DIS kinematics 

Let the incoming and the scattered positron have four-momenta k and k' respec- 
tively. Then the exchanged photon virtuality is 

g2 = -q^ = -{k - k' f = 4EE' cos2(6'j2), (1) 

where 9e is the polar angle of the scattered positron relative to the proton beam 
direction. If the overall ep centre-of-mass energy is s = (fc -I- P)^, where P is the 
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proton four- momentum, then Bjorken x is defined as Q /(2P ■ q). The quantity 

E'l-COsOe 

sx L 2 

represents the virtual photon energy, as a fraction of the positron energy, in the 
proton rest frame. 

The cross sections for unpolarised DIS can be expressed in terms of structure 
functions F2 and Fi. It is usually valid to neglect the contribution from Fi; one 
then obtains: 

in terms of the electromagnetic coupling constant a. The overall structure function 
will be denoted by F2 and the structure function for charm-containing events by 

The 7*p centre-of-mass energy or, equivalently, that of the final-state hadronic 
system, is termed W. 

2. QCD Aspects of Heavy Flavour Production 

Perturbative methods in QCD calculations require the occurrence of a momentum 
transfer that is much greater than the QCD scale parameter Aqcd- We shall first 
outline an approach that is widely used in calculating heavy quark cross sections for 
processes involving a high-p^ scatter. Let a hadron H , containing a heavy quark of 
flavour h, be produced from two initial-state partons i, j. The cross section is then 
expressed as the product of two basic terms: 

<7QCD[i3^hx) f(h^H). (4) 

where <tqcd is the perturbatively calculated QCD cross section, i.e. at the parton 
level. Here the symbol x denotes further parton- level products, and the fragmenta- 
tion factor f{h H) is the probability for h to fragment into H in an outgoing jet. 
The inclusive fragmentation probability for x is taken as unity; if we are interested 
in a particular final-state product X then a further factor f{x^X)\s needed. All 
these quantities are kinematics-dependent. 

In HERA physics, the initial-state partons i, j can be taken as quarks, gluons or 
the incoming photon itself. A parton distribution function (pdf), namely f{p — > j), 
is required for the proton. In resolved photon processes, a pdf /(7 — > must be 
supplied similarly to denote the probability that the photon gives rise to the quark 
or gluon i that interacts. 

Finally, expressions of type are summed over all the relevant parton states. 
Thus for direct inclusive photoproduction of a given D meson one can write 

cr(7p Dx) = ^f{p^ j) cTQCDilj -> cx) f{c ^ D). (5) 



The following points need to be noted: 
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(i) The QCD term (Jqcd may be calculated at LO or to higher orders; 

(ii) The assumption that a process may be split into separate terms in this way 
is known as factorisation. A given fragmentation factor or pdf is normally 
assumed to be independent of the perturbative QCD process ij hx, i.e. it 
is universal. A mention of factorisation usually implies universality. 

(iii) The pdf 's and fragmentation factors incorporate a variety of initial-state and 
final-state processes. For example gluon radiation may occur and, if sufh- 
ciently hard, might legitimately be regarded as constituting a higher-order 
QCD process. By imposing a factorisation scale, one defines up to what mo- 
mentum transfers such a gluon or its products will still be called part of the 
proton structure or part of the jet. 

(iv) In general, these are non-perturbative quantities and must be obtained either 
from models or from experimental measurement. The assumption of univer- 
sality implies that a fragmentation factor measured in one experiment can be 
used correspondingly in another experimental context. The pdf's are usually 
obtained from fits to large collections of experimental data. 

(v) As the order of the QCD calculation rises, the validity of the above assump- 
tions becomes more questionable. 

Behind the kind of scheme outlined above there lies the physical insight that 
hard processes take place over a shorter time-scale than soft processes. Thus the 
production of the heavy quark may be regarded as prior to, and physically distin- 
guishable from, the subsequent slower hadronisation phase. 

2.1. Methods and models 

At this point we present some theoretical models and procedures which are in com- 
mon use in the analysis of heavy quark production. Some further topics will be 
discussed later in connection with the specific areas where they are met. 

• PYTHIA, HERWIG. These well-known Monte Carlo models can calculate a 
large variety of QCD processes at LO. Initial and final state leading-logarithm 
parton showers are incorporated to simulate certain kinds of higher order ef- 
fect. Hadronisation is performed in two different ways. In PYTHIaJ] colour 
strings are constructed between the final-state partons, and then hadronised 
according to a set of phenomenological prescriptions. In HERWIG,! a series 
of parton clusters are produced and allowed to decay. Both these models are 
extensively used, and form the basis of further models. The PYTHIA colour- 
string fragmentation is also known by the name of JETSET. 

• Peterson fragmentation model. A standard representation of the fragmen- 
tation factor f{h H), based on phenomenological considerations, is the 
Peterson formula:! 

fih^H) = PDiz) = P ^ (6) 

z[l- 1/z - e/(l - z)] 
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in which the hadron of interest, H, is produced at a given Feynnian z rela- 
tive to the momentum of the heavy quark h, with z — {E + + p\\)h 
where p\\ denotes the momentum component along the h direction. A is a 
normalisation constant, and P is the total probability for h to fragroCTit to H . 
The so-called Peterson parameter e is determined from experiment ilJ a small 
value means that the fragmentation is peaked near z = 1. For c and b quarks, 
respectively, typical values of e are 0.035 and 0.006. There is a question of 
whether the hadron is produced, on average, along the same direction as the 
quark, as taken in Peterson fragmentation, or whether the colour string from 
the proton remnant might pull the hadron to a higher rapidity (the "beam- 
drag effect"Ll). HI have concluded that the effect is unimportant in DIS charm 
production.lij On the other hand ZEUSlij have obtained an improved descrip- 
tion of charm rapidity distributions when the hadronisation procedure from 
JETSET or HERWIG is used in place of the Peterson method. These issues 
have been discussed further by Bodwin and Harris Jll An analytic fragmenta^ 
tion function with possibly better relativistic properties is that of fowler .Ej 
Others have been given by Kartvelishvili et al.t^ and Collins et al. ;tZl see also 
the review by Frixione et al.lla 

Fixed flavour-number schemes. In heavy quark production up to NLO, two 
main schemes have been proposed. In the fixed flavour-number scheme, often 
referred to in the charm case as the massive charm scheme, the incoming 
photon and proton are given hadronic structures which contain only three 
quark flavours (u, d, s). QCD interactions are then generated which produce 
heavy quark pairs {qhqh), whose dynamics is calculated using a realistic quark 
mass assignment (e.g. nic = 1.5 GeV). The heavy quark then fragments — 
e.g. using the Peterson formula — into an observable hadron. Since heavy 
quark excitation processes (see Section 3) are not treated, and may be impor- 
tant at high energies, these schemes are expected to work best at pt — 0{mh). 

Variable flavour-number schemes. To enable quark excitation to take place, 
charm and beauty are treated as active flavours in the proton or photon. In 
the more common massless ( "zero- mass" ) versions, the mass ruh is treated as 
zero up to the final hadronisation which produces a massive H hadron. QCD 
processes taking place at pT ~ to/i may therefore not be accurately described, 
but this approach should work well at high pT- 

DGLAP evolution. Given a set of pdf's or a structure function which apply 
at low Q^, one may wish to know the corresponding_nuantities at higher 
for a given Bjorken x value. The DGLAP equationglj evaluate the necessary 
gluon radiations and splittings that occur in this transition, using perturba- 
tive QCD to a given order. 

CCFM evolution. As an alternative to DGLAP evolution, the CCFM schemed 
uses a different form of factorisation in which parton densities are explicitly 
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treated as a function of transverse momentum (which is integrated over in the 
DGLAP treatment) and an angular ordering parametei_is introduced. This 
is a development of the earlier BFKL evolution scheme ,Eil which evolves am- 
plitudes over X instead of and claims advantages at small x. With CCFM 
there is a matched transition to the DGcLAP regime; it is implemented in the 
CASCADE Monte Carlo among others-ESy 

• Parton distribution functions. A variety of pdf 's are on the market to describe 
the quark and gluon densities in the proton and the photon. They must be 
provided in versions suited to the order of the QCD calculation to which they 
are to be applied. In the case of the proton, the pdf's are usually obtained 
from fits to large collections of experimental data, and are updated when new 
experimental results are announced. The constraints of DGLAP evolution 
are normally applied. Two groups active in this area areCTEQ and MRST, 
whose publications may be consulted for further details nJ'CJ In addition the 
HI and ZEUS collaborations haxecalculated their own proton pdf fits on the 
basis of the DGLAP equations.E30 

Photon pdf's are normally obtained based on the Vector Meson Dominance 
model (see Section 6), so that the partonic structure of the photon at low 
virtualities uses that taken for virtual mesons. An additional direct coupling of 
the photon to quark pairs ("box diagram") is included, and allowance is made 
for free parameters, such as the gluon density in the photon at low virtuality. 
Evolution in from a base value Qa_is pciformccL The modds cited in 
the analyses quoted here, namely GRV,!^ GSP AFG,eI and SaSP differ in a 
number of technical aspects, such as the value of Qqj the factorisation scheme, 
treatment of heavy quark thresholds, and the number of free parameters that 
have been fitted to published data. The SaS model extends the pd& to virtual 
photons. A review of some of the issues has been given by Vogt.E3 

3. Open Charm Production at HERA 

The production of heavy quarks {qh) in ep collisions takes place through three main 
mechanisms, which are here briefly outlined in turn: 

In photon-gluon fusion (fig. |^a), the incoming virtual photon interacts with a 
gluon from the proton, so as to form an outgoing qnqh pah through the reaction 
7*g — > qhqh. The LO diagram often offers a good description when the quark pair 
emerges at high pT . The fusion process can give a direct measurement of the gluon 
content of the proton; this complements the determinations through QCD fits to the 
proton structure function F2 as a function of . We shall refer to "photon-gluon 
fusion" in the context of photoproduction, and "boson-gluon fusion" in DIS, where 
7* and W/Z exchanges may both be present. 

In excitation processes, the incoming 7* interacts with a heavy quark that is an 
effective part of the proton structure and scatters it out of the proton. This may 
occur through resolved photon processes and through the direct "QCD Compton" 
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Fig. 3. D* T— iD meson mass difference, AM, illustrated from ZEUS photoproduction data 
(preliminary) ,e3 using K-k pairs selected to have invariant masses in the range 1.83-1.90 GeV. 

process, ^*qh ^ whose cross section is somewhat smaller than that of the fusion 
process. The proton has only light valence quarks, and so any c or b content in its 
pdf is generated by DGLAP evolution — specifically, through a gluon fluctuating 
into a qhq_h pair. Similar considerations apply to the photon. 

In production by fragmentation, the heavy quark is found in a jet that originated 
from a light parton: specifically, an outgoing gluon splits into a cc or hh pair. Frag- 
mentation may be treated entirely phenomenologically, obtaining the distributions 
from fits to data, or it may be calculated in perturbative QCD, e.g. in the process 
79 ~* .99 with g qhQh, namely "gluon splitting". This is a higher order QCD 
process than the LO fusion term jg — > gnqh,- 

The distribution of an observed hadron H 'm a, jet requires in any case a phe- 
nomenological treatment of the hadronisation process. 

3.1. Charm meson states 

The lowest charm states are the cd and cu and mesons, which are produced 
along with excited states such as the D* . The D* mesons are easier to identify 
because of their decays to the D states through emission of a pion that is almost 
at rest relative to the D* . The decay chains D*±(2010) ^ D°{imA) + 7r±, with 

— * K^TT^ or D° — > isT^Tr^TT+Tr" (-1- cc.) are commonly selected since all the 
final state particles are charged, and hence accurately measurable. A typical mass- 
difference plot is shown in fig. |[ No peak is seen in the distribution that uses the 
wrong-sign combinations: this quantifies the background to the signal. 

In recent ZEUS photoproduction data, the has been measured both as the 
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Fig. 4. ZEUS data (preliminary) on D meson production from K^tt^ pairs (a) accompanied by 
an extra pion giving a mass difference AM consistent with D* decay, (b) without an extra pion 
of this type. 



product of D* decay and in direct production.cfl Figure ^ shows the Kir mass distri- 
butions for these two cases. After selection on the AM peak, a dramatic Z)° signal 
is seen, while the directly-produced 13° signal lies on a much higher background, 
although the number of events in the peak is much larger. In the latter case the 
wrongly-assigned Kt: combinations remain in the plotted background. From these 
plots, ZEUS have evaluated the ratio Py = V/ {V + PS) representing the production 
of the vector D* (spin-1) relative to the pseudoscalar D. This is predicted to be 
0.75 from simple spin statistics, but the inclusion of decays from heavier charm 
states and the effects of hadronisation can reduce this value. ZEUS obtains the 
value Py — 0.546 ± 0.045 ± 0.028 (preliminary) J23 in agreement with measurements 
at LEP.* 

Other mesons with a single c quark are also observed at HERA. The cs system 
has as its lowest state_the 13^(1969) which has been detected by ZEUS in its 07r^ 
decay mode (fig. ||a).Ea Details of the production characteristics will be discussed 
below. The ratio of the cross sections for and D^* production can give a 
measure of the strangeness-suppression ratio, a phenomenological parameter of the 
Lund string fragmentation scheme. The cross section ratio is found by ZEUS to be 
0.41 ± 0.07^0 05 '^•lO(t'-r-); where the last error is the uncertainty on the (t>n 
branching ratio. This value is in close agreement with the corresponding result from 
e+e^ data. It corresponds to a strangeness-suppression ratio of 0.27 ± 0.04lg Q3 ± 
0.07(b.r.), in agreement with most other experiments and with the standard default 
value of 0.3 in JETSET. Both these and the above results support the concepts of 
factorisation and universality in charm fragmentation. 

A further charm state recently reported by ZEUS is the I3^j^(2536), which decays 

*Here and elsewhere, unless otherwise stated, the first of two quoted uncertainties is statistical 
and the second is systematic. 
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Fig. 5. ZEUS data on cs meson photoproduction: (a) (1969) signal seen in the decay D^. — > 
07r± K+K~-K^ (b) D^j(2536) signal seen in the decay Of^ D*^K° with Z)*± D°tt^ 
and D'^ — > K^ir^ + c.c. (preUminary). The AM between the states before and after the Kg 
emission is added to the standard mass values of the decay products. 



by emission of a K° meson to a This is one of a cluster of L = 1 states of 

the cs system. A clear signal is seen in photoproduction, (fig. ||b). The lyroduction 
rate and helicity characteristics are consistent with other measurements. lZI Evidence 
for the detection of other P-wave charm meson states has also been presented by 
ZEUS.L3 These are neutral states in the mass range 2.4 to 2.5 GeV. Further statistics 
are required to clarify the observations. 

3.2. Inclusive photoproduction of open charm 

ZEUS have published a study in which inclusive differential cross sections for D*^ 
are compared with several "massive" and "massless" charm schemes in NLO 
The KniTTT decay channel of the D provides a check on the higher-quality measure- 
ments using_ftr7r. The massive charm, fixed flavour-number scheme of Frixionc ct 
al. (FMNR)c3 is taken with several values of the Peterson fragmentation pararoiicr 
e, along with the masalfsa, variable flavour-number calculations of Knichl ct aLcJ^H 
and of Cacciari et al.E^uJ The two massless calculations differ in the factorisation 
schemes applied to the fragmentation. 

The shape, although not the normalisation, of the pr distribution of the D*^ is 
fairly well described by all the models, although the description cannot be said to be 
perfect ffig. The issues have been recently discussed by some of the theoretical 
authors.CJ Pronounced differences, however, appear in the r] distributions; these are 
shown in fig. ^ |^ integrated over pT tom different thresholds. In the highest pT 
range^d), the scheme of Kniehl et al.E3 is more successful than that of Cacciari 
et al.,E3 whose predictions are too low. Although the latter scheme does appear to 
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Fig. 6. D* photoproduction: ZEUS distributions in pT compared with models (see text). Outer 
error bars are statistical + systematic combined in quadrature; an overall normalisation uncertainty 
of ±4% is not shown. Note that the data here and in the next figure are quoted as ep cross sections. 



work in the lowest pT range measured (a), this is probably of little significance since 
neither massless calculation is appropriate for pT ~ rric. In fact even the massive 
scheme can be brought to reasonable agreement only by the artificial recourse of 
setting rric as low as 1.2 GeV, while the massless schemes are unconvincing for 
PT < 6 GeV. 

Differences at the 20-25% level can be made by varying the photon pdf, as 
illustrated in fig. |^ - where (e), (g) use the calculation of Kniehl ct al. and (f), 
(h) use that of Cacciari et al. For pT < 6 GeV, the GS-G photon structure gives 
better results for the Kniehl et al. calculation due to its high c quark density, similar 
to that of the u quark. Since the quasi-real photon is believed to couple mainly 
to light-quark hadronic states, this structure may lack plausibility. It cannot, of 
course, be used in the massive-charm, fixed three-flavour-number scheme. 

Overall, the message from these results is that the massive scheme has trouble 
between 3 and 6 GeV, where its chances of working are best, while the massless 
scheme of Kniehl et al. does seem to become more successful as pr increases: higher 
Pt measurements are required to confirm this trend. Results from Hl@ have, also 
been compared with the massive and massless charm calculations of FMNRe3 and 
Kniehl et al.c3 Cross sections as a function of laboratory rapidity are shown in 
fig. ||. At low Pt the results are qualitative similar to those of ZEUS, but at high 
Pt the data lie significantly above the calculations, in strong contrast to those of 
ZEUS. There are apparent differences between the two experiments here which need 
explanation. 

Measurements of Ds productiongivc further comparisons with theory. _[n fig. 
, ZEUS first compare their data,E3 together with the previous D* data,Ej with 
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Fig. 9. Differential cross sections for inclusive Ds production (ZEUS) In the left-hand plot, the 
upper (thinner) curves of each type have rric = 1.2 GeV and the lower (thicker) have rric = 1.5 
GeV. The BKL calculation (right) uses rric = 1-5 GeV. 



the FMNR model at NLO. The fragmentation to the D states was performed with 
a Peterson formula: the /(c — > D*^) and /(c Ds) distributions then differ by 
a simple scale factor. As discussed above, a standard calculation (heavy curves) 
fails to fit the D* data although a low rric value will almost suffice. The Dg data 
are consistent with this, but the statistics are at present pom'. Another type of 
theoretical model has been given by Berezhnoy et al. (BKL),cZl who relate the ex- 
perimental value of /(c D*^) to /(c Dg) by tuning the colour-singlet and 
colour-octet contributions to a description of the hadronisation. The D*^ data are 
used to obtain the relevant singlet/octet ratio, yielding a prediction for the Dg cross 
sections. If both components are used, there is a fair description of the ZEUS data. 

3.3. Charm in dijet photoproduction 

More information can be obtained if the charm meson is associated with the produc- 
tion of jets. ZEUS have presented two analyses on this topic; note that at present, 
the inclusive Z)*-containing events are measured with a defined experimental accep- 
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tance, but with no specific association of the D* with a jet, or information about a 
second charm particle in an event. The two highest-i^T jets found in an event are 
used. In fig. [l0|a, the distribution in x°^^ is compared with predictions obtained 
using HERWIG, the total predictions of the model being normalised to the data. To 
a first approximation, HERWIG describes the shape of the distribution well, despite 
the QCD calculation being basically at LO. Above x°^^ « 0.75 a strong direct peak 
is seen both in the data and the theory. A strong resolved contribution is required. 
The tail to lower x°^^ of the HERWIG direct contribution can be attributed to 
effects of hadronisation and initial-state gluon radiation. The small contribution of 
Z)*'s arising from b quark production and decay is included. 

To obtain this reasonable description of the resolved component, the GRV-G 
HO photon pdf was usedE3 This has given good results in many photoproduction 
analyses at HERA. It is possible to separate between the charm-excitation and the 
other charm contributions (e.g. gg — > cc) in the HERWIG resolved calculation; 
the latter contributions are very small. The charm-excitation here is believed to 
come mainly from the photon: its nature is discussed further in section 3.5. The 
charm excited from the proton ends up mainly at forward angles outside the present 
experiments' acceptanceJ13 

A parton-level comparison with the massive NLO calculation of FMNR is un- 
successful (fig. |lo|b). There is no explicit charm excitation in this calculation, by 
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definition, and three-flavour pdf 's are used. The size of the direct peak is accurately 
accounted for, bearing in mind that gluon radiation pushes some of the events of 
direct origin into the low-x"^^ taiL The effects of hadronisation would further in- 
crease this tail but are estimated to be small. So it is clear that the low-x"^^ data 
are not described by the FMNR model, even with the desperate recourse of putting 
TOc = 1.2 GeV, which would anyway tend to spoil the agreement in the direct peak. 
This suggests as a conclusion that the fixed flavour-number scheme is inadequate 
in the context of charm photoproduction at moderate pr- 

This work has been extended by ZEUS, this time using PYTHIA as the reference 
Monte Carlo, in a study of dijet angular distributions in Z?*-containing events. The 
direct process proceeds by quark exchange, while the resolved processes have a 
mixture of exchanges dominated by that of a gluon. These are respectively spin-i 
and spin-1 exchanges, which have very different distributions in the scattering angle 
0* in the dijct rest frame. Experimentally they can be distinguished, to a reasonable 
approximation, by comparing the direct-dominated, and resolved-dominated event 
samples with > 0.75 and < 0.75 respectively.cS 



Preliminary results in fig. 10; show a steep angular rise in the direct-dominated 
events, in marked contrast to the gentler behaviour of the resolved-dominated 
events. The agreement with expectations is very satisfactory considering that the 
calculations are at LO. A next logical step will be to tag the jet with which the 
charm is associated, and examine the forward and backward distributions of the 
charm-tagged jets in the two classes of event, in order to study charm excitation in 
the photon and in the proton in more detail. 

3.4. Open charm in DIS 

The production of charm in DIS processes is dominated by the LO boson-gluon 
fusion diagram. Both HERA collider experiments have measured the inclusive pro- 
duction of the D*^ in DIS and have used this to measure the part of the proton 
structure function F2 which represents charm production. The method employed is 
to identify the D*^ signal, correct it to an inclusive D*^ production cross section in 
bins of X and , and then use a theoretical model to relate the D*^ cross section 
to F|, defined at the parton level. One writes: 

to obtain an experimental F2 result at bin-averaged values of x, Q^. The theoretical 
terms take account of the heavy quark fragmentation and the finite phase-space 
acceptance of the cross section measurement. The cross sections are taken to be 
sufficiently smooth-varying so that bin migration effects can be neglected. 

HI seek to compare their results with predictions from a DGLAP-based evolution 
scheme and also from a CCFM-based sch£ine.E3 To do this, they first make use of 
the hadronising LO Monte Carlo AR0MAE2I to evaluate inclusive D*^ cross sections 
from the observed event numbers. They then perform the subsequent calculation of 
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Fig. 11. Structure function F,^ (a) measured by HI and ZEUS using an NLO DGLAP scheme to 
correct to full phase space and (b) by HI using a CCFM-based scheme. The preliminary results 
are compared with the respective theories. The shaded band represents the uncertainty on the c 
quark mass (1.3 - 1.5 GeV). Note that the ZEUS and HI data are not at identical values. 



F| in two ways, by means of (a) HVQDIS, which uses DGLAP, and (b) CASCADE, 
which uses CCFM. The HVQDIS programEil uses a fixed flavour-number model for 
DIS in which the boson-gluon fusion process is calculated to NLO. The phase- 
space acceptance correction differs significantly in the two approaches, and so HI 
have presented two different sets of experimental structure functions (fig. |ll|), each 
compared with the absolute predictions of a theory of a type similar to that used 
to perform the correction. The CCFM prediction includes proton gluon densities 
measured by HI in an overall F2 determination using this model. The shaded band 
represents the effect of varying the c mass between 1.3 and 1.5 GeV, which is the 
dominant uncertainty, and in the CCFM case the effect also of using a proton gluon 
density obtained with an NLO fit to F2. 

Overall, the CCFM-bas£d calculation describes the data better than the NLO- 
DGLAP based calculation,E3 the differences appearing at lower values of and x. 
Similar conclusions were obtained when comparing the inclusive D*^ cross sections 
with these two models. 

In a like manner, ZEUS have obtained D*^ cross sections using RAPGAP to 
evaluate the hadronisation corrections and HVQDIS as the parton-level modelEl 
The inclusive cross sections agree well with HI and with HVQDIS, but a JETSET- 
style fragmentation function is required, rather than the Peterson formula, to per- 
form the hadronisation corrections. There are some further technical differences 
from the HI approach. ZEUS also find good agreement between their mea- 
surements and results from an NLO QCD fit using ZEUS proton pdf'sS This is 
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Fig. 12. (a) Charm structure function measured by ZEUS using Et&iNLO DGLAP scheme to correct 
to full phase space, and compared with a fit using this scheme. E3 Uncertainties of around ±10% 
in the normalisation due to the luminosity measurement, the charm meson branching ratios and 
the charm fragmentation function are not included, (b) The ratio of the charm structure function 
to the total structure function i^2- 

illustrated in the F2 results plotted as a function of (fig. |lja) for different values 
of X. The variation of F2 with — namely, scaling violation — is very much in 
evidence, and is well described globally by the theoretical fit. The uncertainty on 
the theory (dashed lines) is dominated by that on the c quark mass. 

The ratio of F| to F2 is presented by ZEUS in fig. [l^b. The fraction of charm 
events is substantial, owing to the flavour-insensitivity of the boson-gluon fusion 
process at high px- The fiattening of the theoretical ratio as x becomes smaller is 
a feature of the boson-gluon fusion model: the rising charm production follows the 
rising gluon density in the proton, accompanying the similarly rising quark density 
which steers i^2- At high x this mechanism should fail, since the gluon density 
falls off and the proton has no charm valence quarks. The experimental data are 
consistent with these theoretical expectations. Overall, therefore, a clear physical 
picture seems to emerge. 

Using a larger event sample from the period 1998-2000, ZEUS have confirmed 
that HVQDIS gives a^good overall description of the D*^ production, up to 
values of 1000 GeV^.EJ An unexpected observation is apparent in the preliminary 
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results, namely that the cross sections measured using an incoming electron beam 
are higher than those using a positron beam, although both sets individually remain 
consistent with HVQDIS. The effect appears highest at high . No explanation 
for the difference is yet on offer, and further investigation is required. 

HI have also extracted the gluon density in the proton on the basis that the 
charm production cross sections are dominated by the boson-glucui fusion process, 
and hence depend strongly on the proton gluon content. Resultsc3 using D* data 
from photoproduction and DIS are plotted in fig. |l^ and are compared with those 
from a QCD fit to F2 data. The F2 measurements depend primarily on the coupling 
of the virtual photon to quarks in the proton, with the DGLAP equations generating 
the gluon distribution in the course of the fitting procedure. The excellent agreement 
seen between the two types of method represents a significant triumph for the 
standard QCD description of the ep system. 



3.5. Photoproduction-DIS relationship 

To probe the nature of the charm process further, ZEUS have compared the pro- 
duction of charm in DIS with that in photoproduction.LJ Dijet events containing a 
D*"^ meson are selected, for « and in the range 1.0 to 5000 GeV^. A value 
of x'^^'^ is calculated for each such event, from which the authors evaluate the cross 
section ratio between events in x°^^ below and above 0.75, i.e. an approximate 
"resolved/direct" ratio. This is plotted as a function of , and compared with 



theoretical predictions in fig. 14. The effects of uncertainties in the proton pdf 



largely cancel in taking this ratio, so that we are mainly looking at the properties of 
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Fig. 14. Ratio of low to high 2;°^^ cross sections of dijet events containing a D* as measured in 
ZEUS (preUminary) . Jet pairs with Et > 7.5, 6.5 GeV are chosen, and D* with pT > 3.0 GeV. 
The results are compared with (a) HERWIG prediction calculated with no charm requirement. 
The possible effect of acceptance corrections is indicated, (b) HERWIG (with charm requirement) 
using the SaSlD photon structure with and without suppression of the photon parton structure 
with Q2_ (b) CASCADE and AROMA, 



the photon. The ratio is found to be approximately constant with Q^, a behaviour 
which contrasts strongly with that found in events lacking a charm requirement, 
and which cannot be accounted for by acceptance effects. Such behaviour is pre- 
dicted by the Monte Carlos HERWIG, CASCADE and AROMA. With HERWIG, 
the SaS-lD photon pdf's were used with and without the option of suppressing 
the hadronic photon contribution with increasing Q^; however, the data are not 
yet able to distinguish between these cases. The AROMA prediction is flat, but 
quantitatively too low, while CASCADE gives a good overall description. 

These observations are consistent with the expectation that, after factorising 
out the proton pdf, the resolved photoproduction of charm takes place mainly 
through perturbative QCD mechanisms as docs the direct process. Thus the "re- 
solved/direct" ratio should not vary rapidly with Q^. This contrasts with the 
non-perturbative, hadronic processes which dominate light-quark resolved photo- 
production and which show a strong fall-off with . In other words, we are talking 
about "anomalous" diagrams here, rather than about a charm-containing hadronic 
photon structure. 

Nevertheless, further statistics are clearly required to distinguish between the 
two SaS models, and data at higher will also help. 
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3.6. Some theoretical comments 

The HVQDIS model gives good results, but neglects the intrinsic charm that can 
evolve in the proton, in other words charm excitation processes, and so could have 
problems for values above the charm mass squared. However the experiments 
do not yet seem to be sufficiently precise to be sensitive to such an effect. Also, 
HVQDIS does not attempt to describe cc formation from the splitting of outgoing 
gluons, and the formation of final-state J/tp mesons has been ignored in these 
measurements; this contribution is not believed to be at an important level. 

Given the merits and drawbacks of the existing fixed and variable flavour- number 
schemes, one might obviously wish for a more advanced type of scheme that com- 
bines the advantages of both. This would require a model which allows the number 
of active flavours in the proton and photon to vary with or Et, while also 
treating the charm quark mass satisfactorily. In a composite approach, each of the 
previous schemes will dominate in its own particular dynamic range, while a suitable 
"matching" is achieved in the intermfidiate region. A LO model of this type was 
provided hjr Aivazis et al. (AC0T).E3 Further developments have been presented 
by CoUinsJ^ who showed that the ideas can be made to work at all orders, and by 
Amundson et al. (ASTW),E3 who report a program of calculations at NLO. The lat- 
ter authors provide a Monte Carlo program which can be used by experimentalists. 
At present ASTW are able to give a good account of the HERA F| measurements. 
Their procedure is less successful with the differential inclusive charm meson dis- 
tributions, owing to a strong sensitivity to an internal scale parameter within the 
model; good qualitative agreement is obtained with the data, but the theoretical 
uncertainties are very large. ASTW state an intention to extend their calculations 
to higher order. 

There are a number of further implementations of these ideas. This approach 
offers a conceptually attractive way forward to a deeper understanding of charm 
production, provided that a sufficiently high order calculation can be achieved. In 
principle it should be applied also to beauty production. 

It has been pointed out by Chuvakin et al.Ej that the fixed flavour-number 
schemes at NLO are very stable as regards scale variations, and are therefore well 
suited for the description of DIS processes over a wide range. On the other 
hand, the cross sections contain large logarithms of the form ln(Q^/m^), and these 
can be resummcd in a variable flavour-number scheme, where the problematic terms 
can be absorbed into thc^sdf 's. Theaa authors demonstrate that variable flavour- 
number schemes (BMSNEI and CSIvS) are able to give good descriptions of D* 
production at HERA. However the results of the calculations are very similar to the 
fixed flavour-number calculations of HVQDIS. In view of the latter's moderate scale 
dependence, there may be at present no advantage in moving to the other schemes, 
at least insofar as the and x description is concerned. 

It should finally be .remarked that considerable care is required in defining 
at the theoretical leveixil For a further discussion of some of the theoretical issues, 
the account by HarrisE3 may be consulted. 
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Fig. 15. (a) ZEUS differential cross section for production of an electron with transverse momen- 
tum p^' relative to the axis of a jet (Ey > 6 GeV) in photoproduction. The calorimeter energy 
uncertainty is shown as a shaded band. The total HERWIG prediction for electron production, 
and for electron production through 5 and b decays, are indicated. Possible backgrounds from DIS 
positrons confined the analysis to electrons, (b) HI measurement of the impact parameter & of 
muons relative to the primary event vertex. Two jets of Et > 5 GeV are required. 



4. Open Beauty 

Both HERA collider experiments have reported measurements of open beauty, i.e. 
the B family of mesons. The latest HI results have benefited from the use of a 
silicon tracking system, and ZEUS will use such a system in future running. In the 
absence of a precise vertex measurement to help identify the B mesons, with their 
relatively long lifetime, the most common identification procedure is by means of 
their leptonic decays. The high mass of the h quark compared with the c means 
that the B meson decay products emerge at higher angles relative to the direction 
of outgoing heavy-fiavour meson or its jet (fig. Important backgrounds come 

from the leptonic decays of D mesons. 

ZEUS have used this method to obtain an inclusive B photoproduction cross 
section, where B denotes all 5-quark_mesons, taking an inclusive branching ratio 
to electrons. HI have useed muons,E3 and have more recently demonstrated the 
effectiveness of a selection on impact parameter (fig. [Tslb), presenting preliminary 
measurements in both photoproductiontJ and DIS.E3 



The results are summarised in fig. 16 1. Since the measurements are defined us- 



ing different experimental acceptances, they are plotted as a ratio of the measured 
cross sections to the predictions at NLO from the HVQDIS program. A Peterson 
fragmentation function was used to calculate B production rates, from which lepton 
cross sections were extracted using the AROMA generator — in other words, a sim- 
ilar method to that used by HI for their D measurements. The HVQDIS program 
yielded a predicted cross section of 11 ± 2 pb for the production of B hadrons at 
HERA, compared with the HI DIS measurement of 39 ± 8 ± 10 pb. AROMA itseff 
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Fig. 16. Open beauty production in various hadronic processes, (a) At HERA using HI and 
ZEUS, as a function of (b) At LEP using 77 collisions with the L3 detector, together with 
open charm results (upper points and curves). Statistical and systematic uncertainties are added 
in quadrature. 



gives a prediction of 9 pb and CASCADE, at LO, gave 15pb. The photoproduction 
theoretical result was calculated using the FMNR program. 

In short, the integrated B cross sections at HERA lie consistently above the 
NLO theory whenever the experimental errors permit a statement. It is therefore 
worth making comparison with results from other types of experiment. At LEP, 
measurements oLcharm and bsauty production have been made in photon-photon 
collisions by LSEZl and OPAL.eSI The acceptance is evaluated using PYTHIA and 
the results compared to NLO predictions from Drees et alB3 The charm situation 
is found to be not dissimilar to that at HERA. L3 jl^b, top) present charm cross 
sections in overall accord with calculations, bearing in mind the c rnass uncertainty. 
OPAL, measuring D* production with pT between 2 Mid 12 GeV,L3 find that the 
massless calculation by Binnewies, Kniehl and KrarneipaO fits the data well, while a 
massive calculation by Frixione, Kramer and LaenenO underestimates the data even 
with rric = 1.2 GeV. The direct process, here referring to 77 qq, is insufficient to 
account for the observed signal, and the authors argue that a resolved photon with 
significant gluon content is required. 

However the beauty rate in L3 is much higher than predicted. The OPAL B cross 
section, selecting on jets with Et > 3 GeV, is in good agreement with L3, whose 



kinematics are similar. In pp co. 
extensive data on B productio: 



7374 



ns at the Tevatron, CDF and DO have presented 
and, again, find that the NLO calculation is an 
underestimate (fig. |l^) .0 However it can be seen that at high pT the disagreement 
is less serious, the largest discrepancies bein found at pT $ 20 GeV, roughly where 
the other experiments have obtained their measurements. This region, then, where 



Heavy Flavour Physics at HERA. 



25 



A 



pp-^bX, Vs=1.8 TeV, ly*l< 1 



CDF 

T+(2S) 




D0 
I Dimuons 

i Muons+Jets 

(Preliminory) 

. Inclusive Muons 



NLO QCD, MRSR2 
Theoretical Uncertainty 



6 7 8 9 10 20 30 40 50 60 70 

p,"" (GeV/c) 

Fig. 17. Beauty production at the Tevatron in CDF and DO, compared with an NLO calculationEI. 
The different points denote different ways of identifying b events; the dashed curves indicate the 
theoretical uncertainty. 



the momentum scale of the hard scatter is of the same magnitude as that of the h 
quark, is where the main trouble seems to lie. 

There are suggestions, however, that a full treatment of the excitation and 
fragmentation contributions within PYTHIA and HERWIG may enable these mod- 
els to fit the Tevatron h data even at lower transverse momenta. It is important to 
generate all possible final states where jet fragmentation may produce B mesons. 
The excitation contributions can be substantial, as indicated in a recent analysis 
from ZEUS (fig. ^8|).ll3 Although the present experimental errors prevent strong 
conclusions from being drawn, a large disagreement with theory is not supported 
by the latter data, although the agreement in the forward direction would be poor 
without the excitation component. 

The value of the h quark mass requires comment. The Particle Data GroupQ 
quote a value centred on 4.25 GeV at production threshold. In_QCD, this mass runs 
with momentum scale, as confirmed by experiments at LEP.LJ However, programs 
such as FMNR are based on a different subtraction scheme which requires the "pole" 
mass, whose value is a little higher, e.g. 4.5 - 5.0 GeV, and is fixed. The Tevatron 
results are accompanied by an NLO curve calculated with rriB — 4.75 GeV^, and 
the dashed curves indicate the uncertainties due to the b mass, the renormalisation 
scale and the Peterson parameter. 

Lipatov et al@ have discussed BFKL-style models which predict the photopro- 
duction cross section for heavy quark-pair production. The essence of these schemes 
is to consider alternative types of gluon distribution to those generated by DGLAP 
evolution. Within the parameter space of this type of approach, these authors are 
able to account for the HI beauty cross section as well as for the total charm cross 
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Fig. 18. Results on open photoproduction of beauty (ZEUS, preliminary) compared with predic- 
tions of PYTHIA with and without the excitation contribution, and with HERWIG. The B signal 
is observed in the muonic semi-leptonic decay mode. The experimental cuts are indicated. 

sections measured by HI, ZEUS and lower energy experiments. A further account is 
given by Baranov, Lipatov and Zotov,Eil who demonstrate that such models are able 
to describe the Tevatron b production cross sections, again provided that suitable 
variants of the theory are selected. 

Although this sounds promising, it is not yet clear that any one implementation 
of the theory is capable of describing the data from all the experiments. To get 
the HI cross section right, the BFKL model requires a &-quark niass value of 
4.25 GeV,E3 which is a little low, while the "semihard" LRSS modelO is good for 
both ZEUS and HI, but overshoots the CDF and DO data if the same parameter 
set is used. The CASCADE Monte Carlo is able to give a good description of the 
Tevatron data,HL3 and, as seen, is fairly satisfactory with ZEUS but not HI. In 
this model, initial-state gluon radiation gives rise to many low-a;°^^ events. The 
authors tend to believe that the HI and ZEUS data are in disagreement. 

At present a summary of this fluid situation is difficult: clearly, 6-quark produc- 
tion provides an testing-ground for a variety of theoretical ideas. As these develop, 
and as the experimental errors from HI and ZEUS decrease, we may expect to 
achieve serious discrimination between the different approaches. 

5. Inclusive Quarkonium Production 

The production of cc and hh mesons has been the subject of extensive study at 
HERA. Specifically, we are referring to the J/ip meson (i.e. the V'(l'S') state), the 
radially excited -0' (or V'(2S')), and the set of analogous T states for the bb system. 
The present section discusses the inclusive production of the J/?/; and '4}{2S) states 
in photoproduction and DIS. Like the production of open charm and beauty, this 
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Fig. 19. Dimuon mass spectra in HI photoproduction.Ej Tlie full 1996-97 data sample is shown 
(top), followed by the selection used for the inelastic J/-(/) measurement and the i/>'(2S) measure- 
ment. The 7p centre-of-mass energy range is 40 - 180 GeV. 



provides an interesting workshop for the study of QCD mechanisms. Experimen- 
tally, the J /ip is conveniently studied at HERA only in its /i^/i^ and e^e~ decays, 
but these are fully adequate to identify the particle with reasonable statistics. 

The other main area of quarkonium physics at HERA is in the diffractive pro- 
duction of vector states qq V through the "elastic" mechanism 7*p Vp. This will 
be discussed in Section 6; at HERA, it is only in these processes that the bb states 
have so far been detected. 



5.1. Production of cc mesons 

In inclusive cc processes at HERA there is no necessity for the bound quark pairs 
to appear in vector states; however only the vector J/tp and ip{2S) states have been 
actually identified, experimentally. Figure ^ illustrates the signals observed in these 
channels by HI. a 

As with the production of open heavy quark systems, these processes are mod- 
elled theoretically in two stages: (a) a hard QCD subprocess, which is calculable 
perturbatively and in which the production of the cc pair is described, and (b) a 
soft hadronisation stage where there is a certain probability for a given hadron to 
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be formed. This may be represented as 

da{-i*p ^ HX) = J2 Ml*V ^ qhqh{n) + x) O" {n) (8) 

n 

in which a series of intermediate qh^h states, labelled n, are considered, and each has 
a certain matrix element (n) for transition into the chosen hadronic state H. The 
perturbatively calculable subprocess acts at a momentum scale given by the mass of 
the qhqh system, or the px of the scatter if higher; it is a short-distance, fast process. 
The intermediate qtqh system may be in a number of angular momentum states, 
and is produced in one of two colour states, namely singlet or octet. If necessary, 
soft gluon radiation ensures the colour-neutrality of the final-state particles. 

A specific implemeiitatiion of this scheme uses the so-called non-relativistic QCD 
or NRQCD approach.E§E3 The initial QCD calculation having first been performed, 
this effective field theory model allows the {n) factors to be expressed in terms 
of known powers of w, the velocity of the heavy quark in the quarkonium system, 
which should not be too large. The matrix elements are otherwise not calculable 
from first principles. However, they should be independent of the original QCD 
process that formed the qq state, and the hope is that by measuring experimental 
cross sections in one context we can obtain information on [n) terms which can 
be applied in another. 

The chief experimental questions at HERA concern: 

• The relative importance of the colour-octet state, or even whether it is im- 
portant at all; 

• The effective value of the mass of the c quark; 

• The need for higher-order QCD calculations. 

There are expectations that the method will work better for hb systems than cc 
systems since the non-relativistic approximations wilLbe more accurate. 

Further details have been given by Frixione et alJ13 and in an extensive review 
by Kramer.EZI This is a situation where several QCD momentum scales apply in 
the same event, and in particular it is assumed that mqV 3> ^qcd- In general, 
the power of v which multiplies the matrix elements depends on the relative size of 
mqV, mqv'^ and h-qcD- 

The simplest version of these ideas is the colour- singlet model, in which colour- 
octet (CO) terms are neglected, since they are suppressed by a power of w'* compared 
to the simplest colour-singlet (CS) term. However they are multiplied by a colour 
factor of 8 and have a different pT dependence, so that their unimportance cannot 
be taken for granted. Theoretically, the octet terms are required in order to remove 
certain infra-red divergences due to soft gluon emission. For best accuracy, one 
would wish the sum in (|^) to include as many terms as possible. In practice, 
it is not possible to determine more than a small number of the (n) factors 
by experiment, and so the hope is that a model with just a few terms included 
will be able to test the assumption of their universality to sufficient accuracy to 
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Fig. 20. (loft) Experimental data from HI and ZEUS on inelastic J/i/) photoproduction, compared 
with theory (Kramer), (right) Data from HI (preliminary) at low z values. 



be usefulH The (different behaviour of the octet from the singlet terms at large 
transverse momentum should also provide a useful handle for our understanding of 
how well these ideas describe the physics. 

Both Hl@ and ZEUsS have presented measurements of inclusive J/V' photopro- 
duction. HI have in addition made the first observations of the -i/;' in this channel. 
Cross sections are given in terms of the of the J/i/'j its centre-of-mass rapidity 
y* ^\ ln(i? —Vz) and its inelasticity z — P ■ P^jP ■ q. Here P, P^ and q 

are the four-momenta of the beam proton, the J/?/; and the incident virtual photon. 
The z parameter is the fraction of the photon energy taken by the J/^ in the proton 
rest frame, and is unity for elastic events. 

Figure ^0|(left) shows the cross sections as a function of 2;, with predictions at LO 
from the NRQCD model. For z > 0.3, the theoretical cross sections are dominated 
by the direct process, with the CO contribution becoming increasingly dominant as 
z increases. The CS resolved cross section is small compared to the CO. The shaded 
band represents the total theoretical uncertainty. The dominant QCD diagram is 
from photon-gluon fusion, — > cc. 

HI have made a special studjH of the low-z production of J/ip^s. Figure 
^(right) shows results that are in excellent agreement with the LO curves, al- 
though the theories quoted here and previously are not identical. The presence of 
the resolved photon contribution seems clearly required. 

It is evident that, while not in disagreement with the LO theory, the data lie 
at the low edge of the predicted range for CS -|- CO, and at the high edge (if a 
similar error band is allowed) for CS alone. The same data are compared to a NLO 
calculation in fig. pl|(a,b), where good agreement is found with a CS prediction 
alone, although the normalisation uncertainties are again large. The theoretical 
uncertainties include those on the c quark mass and the QCD coupling constant a^- 
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Fig. 21. Experimental data from HI and ZEUS on inelastic J/-(/) photoproduction, compared 
with NRQCD theory (Kramer) 
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Fig. 22. HI DIS ep cross sections for various distributions in inclusive J /tj) production (open cir- 
cles), to be compared with the AROMA histogram, and for inelastic {Mx > MVd-GeY) production 
(solid circles), to be compared with NRQCD curves from Fleming and MehenE2l 



have also compared their resuhs with a NLO CS calculationtil and with a 
LO CS calculation supplemented by a LO CO model using fits from CDF and from 
CLEO data. The conclusions are similar: NLO CS succeeds, the others sometimes 
succeed and sometimes fail, but the theoretical uncertainties are large. Clearly, this 
does not yet amount to a powerful test of these ideas, especially for the presence 
of the CO term. The pt distribution (fig. pl]c) demonstrates the importance of 
an NLO calculation. It includes terms which are dominated by t-channel gluon 
exchange and scale as instead of pt^ , but similar t-channel exchanges occur in 
the LO CO contribution. 

In an analysis of DIS J/ip production,^ HI have presented results with and 
without a mass cut which effectively removes most events at z > 0.9, where the 
elastic and CO effects are largest (fi£^|2^). They find that the CO contribution, 
as calculated by Fleming and Mehen,L3 still dominates the remaining phase space. 
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Fragmentation is not included in this LO calculation. The CS term is inadequate to 
describe the data, but including the CO term makes the total too high. HI attribute 
these discrepancies to the matrix elfiinents used in the calculation. Meanwhile 
the Soft Colour Interaction model,o calculated in AROMA, is inadequate at the 
quantitative level, although the shapes of most of the distributions are described 
reasonably well. 

One path to a better test of these ideas is to measure the polarisation of_the 
J/V", where the effects of the CO terms should be more unambiguously visiblcEil It 
would also be interesting to make an approximate separation between the resolved 
and direct contributions, by means of a;^, using the events with a J/ip jet and a 
gluon jet. 

The production of V-''(2S) has been measured by HI relative to that for J/i/? in 
the same kinematic range of inelastic photoproduction.0 The cross section ratio was 
found to be 0.210 ± 0.048 ± 0.032 (preliminary). A substantial fraction of observed 
J/?/' mesons can be direct or indirect decay products of '0'(2S) mesons. Allowing 
for this, HI calculate the ratio of ■(/''(2SI-to directly produced J/ij] as 0.24 ±0.08, in 
agreement with an expectation of O.lS.til 

The NRQCJl model is itself being developed and studied, with a variety of ideas 
on the market.Ej These will in due course invite further experimental test. 

6. Diffractive Processes 

The final topic of this survey concerns the diffractive production of heavy quark 
systems at HERA. As illustrated in fig. such processes occur through the pho- 
ton first fluctuating into a hadronic state V , where V must preserve the quantum 
numbers of the photon. The assumption of vector meson dominance asserts that 
the relevant hadronic states are the members of the vector meson family, namely the 
L ~ 1 qq mesons p, uj, (p, J/i/', f etc. The state V may interact inelastically with 
the proton, but the Optical Theorem requires the existence of an elastic channel, 
in which the fluctuation 7* — > is followed by the elastic process Vp Vp and 
the V emerges on mass shell. A full theory must also take into account excitation 
processes whereby the V or proton may emerge in an excited state. 

From a QCD standpoint, the essential feature of diffractive processes is that 
in the primary scatter, the exchanged object is colour- neutral (colour-singlet). At 
HERA, this condition permits the proton to remain intact, to end up in an excited 
nucleon state, or to become dissociated. However there must be no transfer of colour 
that requires subsequent neutralisation by the generation of a string of hadrons 
along the rapidity range between the forward nucleon system and the remainder 
of the event. This rapidity range usually ends up underpopulated or unpopulated, 
therefore; in fact a so-called "large rapidity gap" is commonly found in the forward 
region of a diffractive event. 

Diffractive processes may be interpreted either in terms of Reggc theory, where 
the exchanged object is a pseudoparticle with the quantum numbers of a vacuum, 
known as a pomeron, or else in terms of QCD models with the typical exchange of 
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Fig. 23. QCD-modelled production of a vector meson in elastic diffraction (from ZEUf: 



two or more gluons (fig. |2^). We shall here, for convenience, refer to the exchanged 
entity generically as a "pomeron" . 

An account of heavy quark production in diffractive ep scattering must place it 
against the backdrop of what is now a broad area of "pomeron" physics. Indeed, 
the manifestations of the pomeron radiated from the proton parallel in many ways 
the behaviour of the photon radiated from the positron. The pomeron is hadronic, 
and it may display parton structure and give rise to "resolved" interactions. In 
the present context, though, we regard it simply as an exchanged object with a 
possibly complex nature, capable of being modelled in various ways. In the elastic 
photoproduction of vector mesons the pomeron exchange is a soft process so long 
as the vector mesons are made of light quarks. However, when cc or bb mesons 
are produced, the process cannot be entirely soft from a QCD viewpoint, and so 
its characteristics change. The situation gives rise to the possibility of perturbative 
QCD (pQCD) calculations. 

In the elastic photoproduction of p, oj and </> mesons, the fluctuation into a vector 
meson is taken to occur before the scatter off the proton. For the vector meson to 
remain intact, the scattering process must then be essentially soft. In the QCD- 
modelled process (fig. |2|), the photon couples to a virtual qq pair, which exchanges 
gluons with the proton, and then finally the emerging qq pair may hadronise into 
a vector meson. In elastic vector-meson production, the coupling to the proton 
selects the latter's gluon content at low mean values of Bjorken x. Since oc 1/a;, 
higher W means lower x, corresponding to a rapidly increasing gluon density in 
the proton pdf (fig. |l3|). One expects to observe in general, therefore, that elastic 
cross sections will increase strongly with W if perturbative QCD is operating. In 
the Regge approach, a weak rise with W is expected. l3 

Thus the important difference between the elastic photoproduction of light and 
heavy vector mesons is that the production mechanism of the former is fundamen- 
tally "soft" at all stages, whereas the latter require a "hard" phase initially. The 
situation changes in DIS. The elastic formation of even a light vector meson from 
an incoming highly virtual photon, with 3> my requires a hard vertex some- 
where, which renders the process in part perturbatively calculable while depressing 
the cross section. 
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6.1. Measurements of vector charmonium 

The interplay of the different hard scales in diffractive vector meson production 
is illustrated in the sels of total elastic cross sections shown in fig. ^ In the 
first two plots, earlierE3 and more recento ZEUS results on p production as^ 
function of W are compared with HI. In the third plot, J/^' data are presented.LJ 
Recent HI p resultai22l are in agreement with ZEUS; the agreement seen in the 
J/ip case is likewise reasonable. The cross sections are presented for the cases of 
photoproduction {Q^ = 0) and for DIS at various values. At a given value of 
the M^-dependence of the total cross sections can be parameterised as , where 6 
is a power to be measured experimentally. 

In the case of p photoproduction, the data over a wide range of W are well fitted 
by the Regge-theory formalism of Donnachie and Landshoff,LZl and indicate S values 
in the range 0.12 to 0.16 at low Q^. As rises, S increases to values in the range 
0.4 to 0.9, but with a large error for the highest data set. The J/tp data sets, 
in contrast, show a large value of 6 already at = 0. A value of S in the range 
0.6 - 0.8 is indicated, with no evidence from HI and ZEUS for variation with Q^. 
The suggestion is that by the time has reached a suitably high value, a common 
type of hard behaviouiunay be found in this W range. 

In the fourth plot|122l J/iP data at = q ar£_combined and compared with 



pQCD predictions using different proton pdf sets,ll£3 all of which give good fits to 
the proton structure function F2 . The data show a clear sensitivity to the pdf set 
— that is, to its gluon component — and suggest a preference for one of them. 
Howeyer this should be treated with caution since the two models chosen, FKS and 



have large uncertainties on their normalisation. They differ in their 
treatment of the gluon distributions in the photon, and in the description of the 
charmonium state. Overall, these models may be said to confirm the shape of the 
variation with W, and the general validity of the pQCD approach to this topic. The 
FKS model describes the data well when renormalised by a factor of « 1.6. 

The coupling of the photon to the vector mesons V is through their quarks; the 
quark content of the p, oj, (f>, J/ip and T mesons predicts total elastic cross sections 
for jp Vp in the respective ratios 1 : ^ : | : | : |, respectively, if flavour 
independence otherwise applies. For five quark types, this is referred to as SU(5) 
symmetry. The uj : p cross section ratio is fo und to be in good agreement with this 
prediction, both in photoproduction and DIsJihI 

However at = both the <j) and J/tp cross sections are lower than predicted. 
This must be attributed to different production mechanisms, even the (j) being 
sufficiently heavier than the p and uj to follow significantly different dynamics. To 
expect SU(5) symmetry to hold, we must operate in a kinematic regime where 
uniform production mechanisms operate, and it would appear necessary that this 
be either one where all the mechanisms are soft, or one where pQCD applies but 
quark mass effects are unimportant. In the latter case this means there must be 
a hard momentum scale considerably greater than Aqcd and the vector meson 
masses. 



Heavy Flavour Physics at HERA. 



36 



T 

^ 0.5 



T 

% 0.; 



T 



f 10- 



(a) 



12 3 4 



-t (GeVT 



ZEUS (prel.) 96-97 
ZEUS 1995 



-t (GeVT 



(b) 



• ZEUS (prel.) DI5 96-99 
■ ZEUS (prel.) BPC 99-00 

* ZEUS (prel.) PHP 
o H1 DIS 95-97 

A ZEUS 95 




[GeVT 



Fig. 25. (a) Ratio of cross sections da/dt for </< and p (upper) and for J/?/) and p (lower) measured 
by ZEUS in proton-dissociative elastic scattering 7p — > VY . (b) Ratio of elastic cross sections for 
J/ij} and p, as a function of photon virtuality as measured by ZEUS and HI. 



A first area to check is scattering at high momentum transfer whe re Z EUS have 
measured proton-(iissociative diffractive production of vector mesons JlH^ The pre- 
dictions appear to work (fig. p5|a): for |t| above approximately 3 GeV or 6 GeV, re- 
spectively, the differential cross section ratios d<T/dt{j*p — > (j)Y) / da/dt{'-f*p — > pY) 
and da/dt{j*p J /if) Y) / da/dt{'-^*p pY) become consistent with asymptotic 
values of 2/9 and 8/9. More high-|t| data are desirable to confirm this to greater 
precision. The corresponding ratio&_£Qr the purely elastic 4> and J/V' production 



have also been measured by ZEUSiliiill they are consistent with expectations but 



the experimental errors are substantial. 

Likewise, the total cross section ratio for relative to p plateaus to a ratio of 2/9 
when exceeds about 6 GeV (which seems to be a hig her v ^li i e than in the case 



of t, indicating that ~t and Q are not equivalent scales) .t2§t23 For J/i^, the data 



are too sparse at present to make a definite statement; however there is a steady 
rising trend with Q ^, a nd an asymptotic value of 8/9 is possible (fig. p5p). 



HI have arguedt22l (fig. 26) that a universal curve can be plotted through all 
the total cross sections for 7*p Vp. The plotted results were taken at a fixed W 
value of 75 GeV and the different cross section values were scaled by the reciprocal 
of their respective SU(5) factors, as listed above. (The T cross section needed to 
be extrapolated in in a possibly questionable way.) After this, a function given 
by Gfit = ai{Q^ + my + 02)""^, with suitable values for the constants a^, passes 
impressively though the data points. 

No explanation was offered by HI for this observation. ZEUS have studied the 
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matter furtherli£il with the conclusion that the postulate of a universal scale variable 
{Q^ +my) is at best an oversimplification. It works for p, to and cj) production, but 
fails for J/ip production when the matter is studied in detail. Taken over a broad 
W range, and at a series of {Q^ + niy) values, the scaled J/ip cross sections are 
consistently higher by around 50% than the light meson data. At least SU(5) is 
broken. 

As well as scattering elastically in diffractive processes, the virtual photon may 
dissociate into an unbound qq system. ZEUS have presented measiixfimmts of 
D* production in such processes, both for photoproduction and DISJli^l'EiHl The 
kinematical distributions, at present with large errors, are approximately r.oTisistont 
with boson-giuon fusion in resol ve d -p om eron predictions using RAPGAPJliil HI 
have measured this process in DISllij and have obtained satisfactory agreement 
with QCD model predictions. In both experiments, a suitable parton model of the 
pomeron is needed, with a dominant gluon component. 

6.2. Regge phenomenology 

The previous sections have treated diffraction principally from a pQCD viewpoint, 
which has an obvious relevance to heavy quark production. Elastic light vector me- 
son photoproduction is described well within the older Regge framework, however, 
and this must now be brought into the discussion. A recent account of this area of 
HERA physics has been given by Abramowicz. 
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Fig. 27. Elastic scattering parameters b and a in the process ■yp — ► J/ipp measured by ZEUS 
(preliminary) and HI. 



The differential cross section for Regge pomeron exchange is given by 

« exp'^'MTTT^ (9) 



d\t\ \WgJ 

■where a{t) ~ a{0) + a't. The constants a, a' are beheved to be universal, and 
characteristic of the pomeron, ■while 60 a-nd Wq depend on the particular process. 
The t distribution depends on 

h = bo+4a'ln{W/Wo). (10) 

A positive value of a' means that the cross section falls more rapidly ■with |i| as W 
increases, a characteristic kno^wn as "shrinkage" . The total elastic cross section, on 
integrating (^), varies as ■where 5 = 4(a(0) — 1 — a' /b). Experimentally one 
finds a(0) k, 1.08 and a' « 0.25 GeV"^ Light vector mesons have 6 w 10 GeV-^. 

There is already a difficulty here in J/il> photoproduction, because a value 6 « 4 
GeV~^ is measured (fig. [2^a), ■which ■would give i5 « 0.1, ■while the measured value of 
5 is around 0.8; in other ■words the actual W dependence is steeper than the Regge 
formula requires. On the other hand, if pQCD ■were the ■whole story, we should not 
expect to see shrinkage in the J differential cross section, since in pQCD the W 
and t dependen ces a. re not coupledJiia From fig. |27| b it is evident that shrinkage is 
indeed observedJl^ This implies non-perturbative physics, ■which is good ne^ws for 
Regge theory. 

Nevertheless, the W dependence represents a problem. In order t o re scue Regge 



theory, Donnachie and Landshoff have proposed a "hard" pomerontiZI as well as 



the familiar "soft" one which describes light vector meson production well.lii3 The 



hard pomeron has q;(0) « 1.4 while the soft one still has q;(0) ~ 1.08. However, 
the experimental t variation gives an intercept of a(0) = 1.193 ± 0.011 ^ooiq- This 
agrees with a BFKL calculation Jlij but with neither the soft nor the hard pomeron 
on its own. 
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Fig. 28. DifTorcntial cross section in \t\ for J/i/' production in DIS with proton dissociation (HI). 
A QCD-based theoretical prediction is shown. 



In a more complex approach one may attempt to combine the two types of 
pomeron. The above authors and HI have shown that a fit involving b oth p o in erons 



is able to describe the observed variation of the cross section with M^.EiiEzj If this 
is becoming a little complicated, perhaps the pomeron is complicated. 

One aspect of simple Regge theory which is found to hold in J /ip production is 
s-channel helicity conservation (SCHC). This has been reported by ZEUS from an 
analysis of the relevant spin density-matrix elements using the angular distributions 
of the J/i/j decay products. §1 

6.3. Proton dissociation 

In diffractive events at high momentum transfers the proton can dissociate 
while the vector meson still emerges elastically, i.e. with no additional associated 
fiua.l-stat p pr oducts Recent ZEUS measurements of this process have already been 



discussed.tHSl In J/ip photoproduction, HI have made a studyliHiJ selecting events in 



which the proton dissociates into a final state with a mass in the approximate range 
1.6-30 GeV, as measured in their forward calormeter system. The value of t is 
obtained from the trans ver se momentum of the measured J /ip. The pQCD-based 



event generator HITVMlHa describes the W and the t distributions well (fig. |2q ). 
The shape of the variation with t appears to be independent of T4^; at given t the 
cross sections rise as with (5 « 1. 

6.4. ip{2S) production 

HI have measured the elastic production of the ipi^S) state both in photoproduc- 
tion and in DIS.E£j'l3 The decay channel ip(2S) — > J/ipn^TT^ was used, and the re- 
sults include the contribution from proton dissociation. A colour-dipole based QCD 



predictionHfl describes the data well (fig. ^9|). There is also a predictionlii3 that for 



m?j/^, the ratio of the elastic tp{2S) and J/^ cross sections should tend to 
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Fig. 29. Differential cross section for tl){2S) production in DIS with proton dissociation: (a) results 
from HI with theoretical comparison, (b) results from ZEUS (preliminary) compared to other 
experiments. 



0.5. ZEUS have measured ip{2S) production through the e+e decay channeLand 
obtain a result that is consistent with HI and with lower energy experimentsJHj 

6.5. Measurements of bb 

The T family are the only bb states measured at HERA, and have been observed 
only in the elastic channel in photoproduction and DIS, th rou gh their decays. 



The dimuon spectrum from HI is illustrated in fig. pO|a.E^ A T enhancement is 
seen with a little difficulty above a background due to the Bethe-Heitler process 
77 —^ fJ-^ IJ-^ , where one of the photons is radiated from the proton. ZEUS have 
observed a similar signal.E£3 The mass resolution is not sufficient to separate the 
T(IS'), (25*) and (35) states, which cover the mass range 9.46 - 10.36 GeV, and so 
the signal-background subtraction is performed over a range that covers all three 
resonances together. 

The cross section is evaluated assuming that the standard branching rate to 
muons applies to all three channels and taking 70% of the signal as T(15'). The 
results from the two experiments are shown in fig. pflb. compared with QCD-based 
calculations by Martin, Ryskin and Teubner (MRT)E23 and by Frankfurt, McDer- 
mott and Strikman (FMS).E^ Within present errors, these models all describe the 
data. A calculation by FKSL23 was found by ZEUS to be too low. 

7. Conclusions 

The HERA collider offers a rich breadth of perspectives through which many areas 
of particle physics can be studied, and in so doing stands as an experimental facility 
unmatched in the world. This is particularly the case in the study of QCD processes, 
where the production of heavy quark systems provides a laboratory for the probing 
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Fig. 30. (a) Dimuon mass distribution in photoproduction from HI. The histogram represents the 
simulated Bethe-Heitler background, (b) T total cross section in photoproduction: HERA results 
compared to theory. 

of many important details of the theory. Mesons containing charm quarks have 
been extensively measured, and a start has been made on the study of beauty 
meson production. These measurements have made a useful start in distinguishing 
between different perturbative and semi-perturbative QCD models. Nevertheless it 
is clear that better experimental statistics at HERA will be exremely valuable in 
the study of charm physics, and indispensible for the study of beauty. The ability 
to study the production of B and T mesons at transverse momenta of the same 
order as the heavy quark mass - yet still "hard" compared with Aqcd ~ can give 
penetrating tests of theoretical ideas. 

Although the available energies are not as high as at the Tevatron, and top 
quark pairs lie outside HERA's reach, the recent improvements to the HERA col- 
lider, together with a variety of enhancements in the experiments, make for exciting 
prospects in this area. Following the machine upgrades, a factor of up to ten in in- 
tegrated luminosity may be achieved by 2005. No serious evidence against QCD has 
been found at HERA, but the phenomenological implementations of the theoretical 
ideas do not yet seem to be perfect. An important onus lies on the theoreticians to 
provide an understanding of these areas of QCD that will match the quantity and 
quality of the promised experimental data. 
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